A high-accuracy optical clock via three-level coherence in neutral bosonic 88 Sr 
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An optical atomic clock scheme is developed that utilizes two lasers to establish coherent coupling 
between the 5s 2 x 5o ground state of 88 Sr and the first excited state, 5s5p 3 Po- The coupling is 
mediated by the broad 5s5p 1 Pi state, exploiting the phenomenon of electromagnetically induced 
transparency. The effective linewidth of the clock transition can be chosen at will by adjusting the 
laser intensity. By trapping the 88 Sr atoms in an optical lattice, long interaction times with the 
two lasers are ensured; Doppler and recoil effects are eliminated. Based on a careful analysis of 
systematic errors, a clock accuracy of better than 2 x f 0~ 17 is expected. 

PACS numbers: 42.62.Eh, 32.80.-t, 32.70.Jz, 42.62.Fi 



Armed with superior resonance quality factors, optical 
atomic clocks based on single trapped ions or a collection 
of laser cooled neutral atoms are expected to outperform 
microwave-based atomic clocks in the near future |J,l2l,b|- 
While a large number of quantum absorbers in a neutral 
atom system provide an advantage in the enhanced short- 
term frequency stability, atomic motion during the probe 
phase seriously limits the attainable accuracy. A compet- 
itive proposal is to localize neutral cold atoms spatially 
in a Lamb-Dicke regime while the trapping potential is 
designed such that its presence does not influence the 
clock transition frequency Q . This scheme can be poten- 
tially realized using a far-off-resonance dipole trap oper- 
ating at a wavelength where the ground and the excited 
state of the clock transition experience exactly the same 
AC Stark shifts. In particular, if the atoms are confined 
near the ground vibrational levels in an optical lattice, 
frequency shifts associated with effects such as Doppler, 
collision, and recoil will all be reduced to negligible levels. 

To realize such a scheme, it is of essential impor- 
tance that the polarizabilities of the two clock states are 
matched to a high degree of accuracy at the "magic" 
wavelength [j| . States with scalar polarizabilities are pre- 
ferred, to avoid the problem of the complex and some- 
times uncontrolled light polarization inside an optical 
lattice. The fermionic strontium isotope, 87 Sr, offers a 
nearly satisfactory solution 0, . The nuclear magnetic 
dipole moment of 87 Sr makes the 5s 2 1 So - 5s5p 3 Po tran- 
sition weakly dipole allowed, with a predicted linewidth 
of about 1 mHz 0, HI • The two clock states have very 
small electronic angular momenta (due to hypcrfine mix- 
ing) and their polarizabilities are nearly scalar. However, 
the hyperfine structure can cause a clock frequency shift 
(possibly as large as 10 Hz) by fluctuations of the lattice 
light polarization. The large nuclear spin (I = 9/2) also 
brings complexity in state preparation and field control. 

In this Letter, we propose a scheme for an optical clock 



that is based on states of true scalar nature, namely 
5s 2 - 5s5p 3 P Q of 88 Sr (with 1 = 0). A direct tran- 
sition between these two states is of course completely 
forbidden, and the clock scheme is based instead on three- 
level quantum coherence established by two probe laser 
frequencies. We note that recently a scheme has emerged 
based on four levels and three probe lasers 0. An esti- 
mate of the expected accuracy for this scheme is not yet 
available. A scheme similar to what is discussed here was 
investigated for trapped ions . 




FIG. 1: Three-level coherence in neutral 88 Sr. 5s 2 1 So and 
are the two clock levels. By applying the probe and 
the dressing laser as shown, a coherent coupling is established 
between the clock levels, mediated by the 5s5p 1 Pi level. 

Our proposal rests on the three-level lambda-type 
scheme shown in FigJIl Particularly promising in this 
context is strontium 0, 0, 0, llij , which we consider here 
in detail. A similar scheme is conceivable for bosonic 
isotopes of other alkaline-earth species (including ytter- 
bium). Since the ground state (5s 2 1 Sq) — subsequently 
denoted by the state vector |1) — and the first excited 
state (5s5p 3 Po) — \2) — of 88 Sr are characterized by a to- 
tal angular momentum of J = 0, electric and magnetic 
one-photon transition matrix elements between these two 
states vanish to any multipole order. \2) decays radia- 
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tively via E1M1 two-photon emission after a lifetime of 
a few thousand years 0. For practical purposes, it is 
therefore legitimate to set the decay width of |2) equal 
to zero. The 5s5p *Pi state, which is referred to as |3) in 
our scheme, can be reached from |1) via El one-photon 
absorption using a probe laser at a wavelength of 461 nm 
(frequency u) p ). One Ml photon from a dressing laser at 
1354 nm (frequency u>d) is needed to resonantly couple 
|2) and |3). The latter decays primarily to |1). The decay 
width is 7 = 2n x 32 MHz. 

Assuming linear polarization of both the probe and the 
dressing lasers, and assuming that the electric component 
of the probe laser (amplitude £q) and the magnetic com- 
ponent of the dressing laser (amplitude Bo) are aligned 
along the same axis, level 3 is characterized by m = 
for our purposes. Let Q p = (3| Dq |1) £q/H denote the 
Rabi frequency associated with electric-dipole coupling 
between levels 1 and 3. Similarly, fl d — (3 \ M |2) Bo/h is 
the Rabi frequency referring to the dressing process. Q p 
and Qd may be taken to be real. 

We first analyze the dressed-atom problem within the 
subspace spanned by |2) and |3) in the absence of the 
probe laser. The equations of motion in the interaction 
picture read 



ic 2 (t) 



d JS d t 



Ca(t) , 



ic 3 (*) = ^e- i ^c 2 (t)-ijc 3 (/). 
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where C2(t) and c^{t) represent probability amplitudes 
and the rotating-wave approximation has been applied. 
The dressing laser detuning is 6 a — uj d — (uj 3 — oj 2 ) [u>i 
stands for the unperturbed frequency of \i), i = 1,2,3]. 

Due to the action of f^, some |3) character is admixed 
to 1 2} (and vice versa), so that a long-lived dressed state, 
1 2), emerges whose decay rate can be adjusted through 
the dressing laser intensity. It is easily seen that to lead- 
ing order with respect to 0,4/(64 + 17/2) 



|2> = {|2>, 



|3> e 



-i(u 3 +Sd)t 
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Here. 
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It follows from the equation for /3 that for small detuning 
the decay rate of |2) is 7 = f^/7. The dressed state 



deriving primarily from |3) decays with a rate that is 
virtually identical to 7. 

Using the atomic structure code described in Ref. 0, 
which accurately treats electron correlation and spin- 
orbit effects in the valence shell, we obtain for the mag- 
netic coupling matrix element |(3|A/o|2)| = 0.022 /j,b- 
Hence, if an effective width 7 of, say, 1 mHz is desired 
(similar to the 87 Sr case), a dressing laser intensity of 
3.9 mW/cm 2 must be chosen and is used for numerical 
results presented in this paper. (By decreasing the inten- 
sity, the effective linewidth can be made even narrower.) 
Under these circumstances, the parameter b in Eqs. (0) 
and (@J, which describes the degree of admixture of |3) 
character in the long-lived dressed state, has a magnitude 
of less than 10~ 5 . In other words, the dressed states in 
this scheme are hardly distinguishable from bare atomic 
states. We will exploit this fact later when estimating 
level shifts caused by sources not included in the simple 
three-level picture of Fig. 

Since there is no direct coupling between |1) and |2), 
it is apparent from Eq. J3J that the transition associated 
with the long-lived dressed state can only be probed from 
the ground state by tuning lo p in the vicinity of UJ3 + 64 — 
oj\, i.e. in the vicinity of the broad absorption profile of 
level 3. Hence, 6 P = u> p — (u>3 — u>i) ~ 84. The long- 
lived and the short-lived states are excited in a coherent 
fashion by the probe laser. They interfere and give rise 
to a narrow dip structure in the absorption profile — an 
effect known as electromagnetically induced transparency 
(EIT) 12]. This dip is more than 10 10 times narrower 
(7 = 2tt x 1 mHz) than the bare absorption profile of 
level 3. The broad background is therefore essentially 
constant across the narrow structure. 

This quantum interference may be treated, to a first 
approximation, within a wave function-based approach, 
but it fails to take into consideration the repopulation 
of 1 1) via spontaneous emission from |3). More impor- 
tantly, the loss of coherence between levels 1 and 2 cannot 
be described. Decoherence is caused by lattice-induced 
spontaneous emission, atomic tunneling between adja- 
cent sites in the optical lattice and ensuing collisions be- 
tween the atoms. Another factor is the loss of phase co- 
herence between the probe and the dressing laser. Both 
aspects must be included in a realistic description. Utiliz- 
ing standard density-matrix theory |l3j . the absorption 
rate from the ground state is found to be 
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where A = (u> p — u>d) — (1^2 — <^i), and j c stands for the loss rate of coherence between |1) and |2). 



If the dressing laser is off, or if -f c 3> 7, the simple 
Lorentzian profile connected with the transition from |1) 
to 1 3) is recovered. Otherwise, if we disregard 7 C for a 
moment, the absorption rate vanishes at A = 0. Note 
that the occurrence of a zero in this absorption profile is 
clear from the connection of this problem to the Fano 01 
lineshape for one bound state (level 2 upshifted by the 
energy of one dressing laser photon) embedded in one 
"continuum" (the 32 MHz wide |3) state). A deep, nar- 
row dip occurs in the absorption profile [lfij . and the 
width of the EIT dip for 7 C = is 



A FWH M=7{l-8^+o(|)} 



(6) 



This demonstrates that not only the position of the EIT 
dip, but also its width is very insensitive with respect to 
fluctuations of 84. A variation of the dressing laser fre- 
quency by as much as 100 kHz induces a relative change 
of the EIT width of less than 10" 7 . The effect is illus- 
trated in Fig. |2(a) (7c = 0). Hence, by measuring W a bs 
as a function of oj p — tod, the clock frequency ll>2 — ^1 can 
be determined with high accuracy by tracking the posi- 
tion of the dip. Neither tv p nor lu^ must be particularly 
stable, only the difference needs to be stabilized. 

Also shown in Fig. [2(a) is the absorption profile near 
A = for relatively large detuning (5d — 5 x 7). Over- 
all, the signal amplitude is drastically reduced in this 
case. There is a narrow absorption peak in the immedi- 
ate vicinity of A = 0, but, as can be seen in Fig. [2(b), this 
peak all but disappears as soon as decoherence becomes 
appreciable. (We have chosen j c = 7 in Fig. 12(b).) On 
the other hand, for small detuning the EIT dip is still 
clearly visible, even though it is not as deep as in the 
absence of decoherence and it has become broader. The 
position of the EIT minimum and the insensitivity of the 
EIT width with respect to laser frequency (5d or S p ) fluc- 
tuations of 100 kHz or so are essentially unaffected. 

The discussion has so far covered only resonant cou- 
plings within the span of |1), |2), and |3). In order to 
evaluate the potential accuracy of the proposed clock 
scheme, we need to consider the AC Stark shifts due to 
nonresonant electric dipole coupling of |1), |2), and |3) to 
states inside and outside the three-level subspace. Our 
calculation of these shifts employs accurate experimental 
and theoretical data on excitation energies and oscillator 
strengths of atomic Sr 0,0,0]. For a probe laser inten- 
sity of 10 fiW/cm 2 , levels 1 and 2 are AC Stark-shifted 
by less than 1 mHz. The dressing laser, operating at an 
intensity of 3.9 mW/cm 2 , causes an AC Stark shift of lev- 
els 1 and 2 of -41 mHz and -20 mHz, respectively. Thus, 
the clock frequency is systematically shifted by about 
+19 mHz. Since the intensity of the dressing laser can 
be stabilized to better than 1 %, we conclude that the 
AC Stark effect due to the probe and dressing lasers can 
be experimentally characterized at the sub-mHz level. 
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FIG. 2: Absorption rate from the ground state of 88 Sr in 
the presence of the two laser frequencies indicated in Fig. 
The laser-induced width, 7, is 2n x 1 mHz. a) Loss of phase 
coherence between levels 1 and 2 is neglected, 7 C = 0. Note 
that the curves for Sd = and <5d = 0.01x7 = 27rx 320 kHz 
are virtually indistinguishable, b) Macroscopic decoherence 
and laser-induced decay are comparable, 7 C = 7. 



To allow for a sufficiently long interrogation time and 
eliminate systematic frequency shifts associated with 
atomic motion, the atoms will be confined in a Lamb- 
Dicke regime in an optical lattice at the magic wave- 
length for the clock transition, as shown in Fig. 3(a). 
The lattice trapping field causes, to leading order, no net 
frequency shift of the clock transition as the level shifts of 
1 and 2 (150 kHz well depth in Fig. 3(a), under a lattice 
laser intensity of 10 kW/cm 2 ) are exactly matched. (The 
studies in Ref. [5j suggest that higher-order AC Stark ef- 
fects are negligible at this lattice laser intensity.) Level 
3 in our scheme, however, is shifted relative to levels 1 
and 2, and this shift depends on the local intensity an 
atom experiences in the lattice. We calculate a relative 
shift of -250 kHz (maximum in magnitude) , at a lattice 
laser intensity of 10 kW/cm 2 and a magic wavelength of 
813.5 nm 0. Figure 3(a) shows the lattice-induced shift 
of level 3 in comparison to 7. This shift is similar — in 
effect and in magnitude — to the drift of the probe and 
dressing laser frequencies. It has, therefore, an equally 
negligible effect on the properties of the EIT dip. 

To make the final estimate of the clock accuracy, the 
level shifts shown in Fig. 3(a) contain the AC Stark shifts 
due to the simultaneous presence of oj p and u>d laser fields 
as discussed above. The lattice field induces a sponta- 
neous emission limited linewidth of the clock transition 
at ^15 mHz. The detuning and decoherence- (laser fre- 
quency variations, atomic motion inside lattice, sponta- 
neous emission, etc)-related changes in the EIT linewidth 
may lead to a small shift in the experimental search of 
the linecenter, but the shift magnitude is estimated to 
be under 10~ 18 . The blackbody radiation induced clock 
frequency shift is estimated to be -0.54 Hz at 300 K, 
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bilization, implemented before ui p and tod interact with 
the atoms. u> p and u>d may fluctuate on the order of 100 
kHz, but their difference fluctuates less than 1 Hz. 

The frequency of LOi oca i can be precisely scanned with 
reference to a stable cavity mode, which of course has an 
inevitable slow frequency drift due to the material nature 
of the reference. However, precise scanning of lo p — luj 
through the EIT resonance enables one to activate a slow 
feedback loop to stabilize lu p — ojd to the value of u>2 — uji , 
i.e., the quantum resonance corrects the long-term drift 
of the cavity resonance according to the clock transition 
frequency. The feedback is applied directly on the 698 
nm local oscillator (uJiocai) to ensure A = 0. 

The tunability of the linewidth, the exquisite insen- 
sitivity with respect to light polarization in the optical 
lattice, the straightforward state control, and the small, 
controllable AC Stark shifts promise that our EIT-based 
clock scheme using bosonic 88 Sr will represent a practi- 
cal and robust approach for optical atomic clocks with a 
competitive performance in both stability and accuracy. 
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Yelin, and M. Lukin for stimulating discussions. Funding 
is provided by NSF, NIST, ONR, and NASA. J. Ye's 
email Ye@jila.colorado.edu. 



FIG. 3: (a) Lattice-confined atoms for use in the EIT clock 
scheme. The AC Stark shift due to the dressing laser on the 
clock transition is shown, as are the level shifts associated 
with the lattice trapping field. The expanded view of the 
Pi state compares the lattice-induced level shift versus the 
linewidth 7. Atoms are cooled to the Lamb-Dicke regime, pre- 
dominantly occupying the ground vibrational level in lattice, 
(b) Experimental implementation of the EIT-based optical 
clock. A pre-stabilized narrow-linewidth local oscillator at 
698 nm establishes a tight frequency track between the probe 
and dressing lasers. This enables a probe of the narrow EIT 
resonance, which in turn provides a long-term clock reference 
for the 698 nm local oscillator. 

with a temperature coefficient of -7.5 mHz/K. Finally, 
the strict scalar nature of the clock transition reduces 
the light polarization and stray magnetic field-induced 
frequency shift below 10~ 18 . In view of the available 
control accuracy of the laboratory temperature and the 
dressing laser intensity, the overall accuracy limit is esti- 
mated to be smaller than 8 mHz, or 2 x 10 -17 . 

Experimental implementation of this EIT-based clock 
scheme is shown in Fig. 3(b). A narrow linewidth (< 1 
Hz) laser at 698 nm (wiocai) will be developed by presta- 
bilizing LOiocai to a stable, passive optical cavity. The 
difference frequency, uj p — u>d, between the probe laser 
and the dressing laser, generated by a nonlinear optical 
crystal, will be phase locked to the 698 nm local oscilla- 
tor. One of the two lasers, such as uJd, as shown in Fig. 
3(b), can be used as the slave to accomplish the optical 
phase lock loop. Thus, u> p — UJd = ^locai- This is presta- 
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